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Abstract 

Background and objective: Metastasis is one of the important reasons for the poor prognosis of hepatocellular carcinoma 
(HCC), abnormal glycosylation plays a pivotal role in HCC metastasis. The goal of this study was to screen and validate the 
transcriptional profiling of glycogenes associated with HCC metastasis. 

Methodology: The differentially transcribed glycogenes were screened out by the Human Glycosylation RT 2 Profiler PCR 
Array, and were identified by qRT-PCR in human HCC cell lines and their orthotopic xenograft tumors. Further analyses were 
performed with K-mean clustering, Gene Ontology (GO) and ingenuity pathways analysis (IPA). Four differentially 
transcribed glycogenes were validated in clinical cancer specimens by qRT-PCR. 

Results: A total of thirty-three differentially transcribed glycogenes were obtained by comparison the transcription in the 
metastatic human HCC cell lines (MHCC97L, MHCC97H and HCCLM3) with the transcription in the non-metastatic HCC cell 
line Hep3B. Seven differentially transcribed glycogenes were selected to further identification in human HCC cell lines and 
their orthotopic xenograft tumors. According to their trends by /(-mean clustering, all of the differentially transcribed 
glycogenes were classified in six clusters. GO analysis of the differentially transcribed glycogenes described them in 
biological process, subcellular location and molecular function. Furthermore, the partial regulatory network of the 
differentially transcribed glycogenes was acquired through the IPA. The transcription levels of galnt3, gcnt3, mania! , 
mgat5b in non-metastatic and metastatic HCC clinical cancer specimens showed the same changing trends with the results 
in human HCC cell lines and their orthotopic xenograft tumors, and the divergent transcription levels of gcnt3 and mgat5b 
were statistically significant. 

Conclusions: The transcriptional profiling of glycogenes associated with HCC metastasis was obtained and validated in this 
study and it might provide novel drug targets and potential biological markers for HCC metastasis. 
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Introduction 

Hepatocellular carcinoma (HCC), which is the sixth common 
neoplasm, ranks the third in cancer mortality in the world [1]. 
80% of the HCC are mainly found in eastern Asia and sub- 
Saharan Africa [2]. HCC is a complex process mediated by 
multiple genes. The risk factors for the development of HCC, 
which interact and cooperate with each other, could increase the 
probability of HCC tumorigenesis. Most people suffering from 
HCC will die within one year after its detection [3]. One reason 
for the high mortality can in part be attributed to extrahepatic 
metastasis [4] and the bottleneck in treatment of HCC is to 
prevent extrahepatic metastasis. The understanding of the gene 
transcription profiling underlying HCC metastasis will provide us 
new theoretical basis for HCC diagnosis and treatment. 



Glycosylation, which can be found in a variety of physiological 
and pathological events, is one of the most important kinds of 
posttranslational modifications. More than 50% of proteins in 
nature are presumed to have undergone glycosylation [5]. These 
glycans not only alter the structures and functions of glycoproteins, 
but also are crucial for cell adhesion and cellular signal 
transduction. Other than that, aberrant glycosylation also plays 
a key role in the underlying mechanism of a variety of diseases [6]. 
Glycans are formed by the catalytic activity of enzymes such as 
glycosyltransferases and glycosidases. The alterations in transcrip- 
tion and translation levels of enzymes are related to corresponding 
changes in the glycan branched structures. Currendy, studies of 
the aberrant glycosylation in HCC have been paid much attention 
and certain achievements have been made. mgat3 and mgal5 were 
proved to be over-expressed in human HCC cell lines, knockdown 
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of mgat3 could promote cells invasion and increase the resistance 
to 5-fluorouracil in vitro [7]. However, the silence oimgal5 in cells 
could inhibit invasion and increase sensitivity to 5-fluorouracil 
in vitro. mgat4a was vital to tumor migration and metastasis 
through altering the giycosylation of CD 147 [8]. It was also 
reported that transcription levels of mgat4b, and mgat5 were 
increased in mice with HCC [9]. Recently, further analysis 
revealed that mgat5 could partially decrease cell adhesion and 
promote cell proliferation through RPTPk [10]. 

In this study, we obtained and identified the transcriptional 
profiling of glycogenes in human HCC cell lines with or without 
metastasis potential and orthotopic xenograft tumors by PCR 
Array and qRT-PCR. The differentially transcribed glycogenes 
were classified and described by K-mt&xi clustering and GO 
analysis, respectively. The possible regulatory network of them was 
built by the IPA. Furthermore, the differendy transcribed 
glycogenges were validated in clinical cancer specimens by PCR 
Array. This work might provide the novel potential drug targets 
and the potential biological markers of HCC. 

Materials and Methods 

Ethics Statement 

Access to human tissues complied with both Chinese laws and 
the guidelines of the Ethics Committee and this study was 
approved by the Research Ethics Committee of Zhongshan 
Hospital and First Affiliated Hospital of Guangxi Medical 
University. All participants have given written informed consent. 

Cell Culture 

The Human HCC cell lines with different metastatic potentials 
(MHCC97L, MHCC97H and HCCLM3) have the same genetic 
background and were established in our Liver Cancer Institute 
[1 1—14]. All of them were cultured in Dulbecco's modified Eagle's 
(DMEM) medium (Gibco, USA) supplemented with 10% fetal 



Table 1. Sequences of qRT-PCR primers. 





Glycogene 


Primer Sequence 


clgaltl 


F: 


5'-ATGAGTGGAGGAGCAGGATA-3' 




R 


5'-TCTGGCACAAAGGGATGAA-3' 


galnt3 


F: 


5 '-CTGCCTCTCCAGGCAACG-3 ' 




R. 


5'-GTAGTACCTGGCGGGTGGG-3' 


gcnt3 


F: 


5 '-AGATGTTG AATGGG AGG AATAG-3 ' 




R. 


5'-TGTTGGTTAGGTGTAATGTGTCTC-3' 


mania! 


F: 


5 '-CTCTCAGCCTACTATCTGTCTGG-3 ' 




R 


5 ' -C AGTTCCTTCC AATACC ACTTT-3 ' 


mgat4a 


F: 


5 '-GCAACGGA AGAACAGGAGTTTC-3 ' 




R 


5 '-AGGTTGGCTACAACACCATGT-3 ' 


mgat5 


F: 


5 '-GCTGCCCCTGTAGG AGAC-3 ' 




R 


5'-GAATCAAGGACTCGGAGCAT-3' 


mgatSb 


F: 


5 '-ATG AGTGGAGG AGCAGG ATA-3 ' 




R 


5 '-TCTGGCACAAAGGGATG AA-3 ' 


st3gal1 


F: 


5 '-CACGAATGGCGTTGGTCTAC-3 ' 




R 


5'-CTCAATCAAAAGGGATGGCA-3' 


fi-Actin 


F: 


5 '-CATGTACGTTGCTATCCAGGC-3 ' 




R 


5 '-CTCCTTAATGTCACGCACGAT-3 ' 



doi:1 0.1 371 /journal.pone.01 07941 .t001 
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bovine serum (FBS, Gibco, USA). Hep3B, the non-metastatic 
HCC cell line (ATCC number HB-8065), was generously 
provided by Cornell University and was cultured in minimum 
essential medium (MEM) medium (Gibco, USA) supplemented 
with 10% FBS (Gibco, USA). All of these cells were incubated at 
37°C with a 5% CO2 in air atmosphere. 

HCC orthotopic xenograft tumor model in nude mice 

The male BALB/C nude mice (5 to 6-week-old) we obtained 
were from Shanghai Institute of Materia Medica (Chinese 
Academy of Sciences, Shanghai, China). The orthotopic xenograft 
tumor models were established as described previously [15-17]. 
Simply, lxlO 7 Hep3B, MHCC97L, MHCC97H or HCCLM3 
cell lines were injected subcutaneously into the upper left flank 
region of nude mice, respectively. After tumor formation, the 
tumors were cut into 2 mmx2 mmx2 mm sized pieces. Nude 
mice that were anaesthetized with pentobarbital 45 mg/kg by 
intraperitoneal injection after regular disinfection were dissected 
and then the tumor tissue lumps were implanted into livers of nude 
mice. The mice were slaughtered at the 35th day after tumor 
implantation and the orthotopic xenograft tumors were immedi- 
ately placed in liquid nitrogen. All of the in vivo experiments 
carried out were strictly complied with the protocol which was 
approved by the Shanghai Medical Experimental Animal Care 
Committee (Permit Number: 2009-0082). 

Clinical cancer specimens 

Fifteen HCC patients under-going resection in 2012 in 
Zhongshan Hospital, Fudan University (Shanghai, China) and 
First Affiliated Hospital of Guangxi Medical University (Nanning, 
China) were enrolled in this study, including ten patients with non- 
metastatic HCC and five patients with metastatic HCC. All the 
tissue samples collected from the patients were performed to qRT- 
PCR analysis. The detailed information of these fifteen patients 
was described in Table SI. 
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Figure 1. Hierarchical clustering. (A) Hierarchical clustering of 84 glycogenes. Each row represented a single gene; each column represented a 
cell line. The expression levels of glycogenes, which were the average value of three test results, were shown by the color scale: red represents a 
target gene with high Ct value while green represents a target gene with low Ct value. (B) Hierarchical clustering of 3 glycogenes acted on 
glycosphingolipids. (C) Hierarchical clustering of 16 glycogenes generated and altered mature N-linked glycans. (D) Hierarchical clustering of 18 
glycogenes generated and altered mature O-linked glycans. 
doi:1 0.1 371 /journal.pone.01 07941 .g001 



Human Glycosylation PCR Array and qRT-PCR 

Total RNA was harvested from cultured cells, orthotopic 
xenograft tumors and clinical cancer specimens with TRIzol 
(Invitrogen) and purified with the RNeasy MinElute Cleanup Kit 
(Qiagen) according to the manufacturer's instruction. After 
quantification using a Nanodrop spectrophotometer, 10 |lg total 
RNA was reversed transcribed into cDNA using the Revert Aid 
First Strand cDNA Synthesis Kits (Fermentas), respectively. 
Besides, cDNA prepared respectively with Super Array PCR 
master mix (Cat. No. PA- 11 2) was needed in order to perform 
PCR Array. 

The Human Glycosylation RT 2 Profiler PCR Array, including 
84 key genes encoding glycan processing enzymes i.e. glycosyl- 
transferase and glycosidase for several important sugars: galactose, 
glucose, mannose, N-acetylgalactosamine, N-acetylglucosamine, 
fucose and sialic acid, was designed and produced by the Qiagen 
company (German). 

The PCR Array and qRT-PCR were performed on IQ5 
machine (Bio-Rad) and the PCR cycling condition was set as 
follows: 95°C for 5 min, 40 cycles of 95°C for 15 s, 60°C for 15 s 
and 72°C for 20 s. Each test was run three times and the mean 
was taken to eradicate any discrepancies. 

The used primer sequences in qRT-PCR were summarized in 
Table 1. 

Statistical analysis and database searching 

Take advantage of the Ct obtained from PCR-Array, the MeV 
4.8.1 was used to map the hierarchical clustering of all the 
glycogenes in the PCR Array, the differentially transcribed 
glycogenes generated and altered N-linked glycans, O-linked 
glycans and glycosphingolipids, respectively. 

The fold change was normalized to the expression of the 
housekeeping gene, fi-Actin, and was calculated respectively for 
each cell using 2~ AACt , where AACt = (Ct target gene- fi-Actin) 
MHCC97L/MHCC97H/HCCLM3-(Ct target gene-fi-Actin) 
Hep3B. Glycogene with differences greater than 2-fold (PS0.05) 
was considered as the differentially transcribed glycogenes and 
they were divided into six clusters by a ^-means approach, which 
is a kind of clustering algorithm for grouping genes or proteins 
with a similar expression pattern [17,18] and the GO was used to 
describe. IPA (Ingenuity System Inc, USA) is a proof of knowledge 
based comprehensive software of data analysis that can help 
researchers model, analyze, and understand the complex biolog- 
ical and chemical systems in life science research [19,20]. The 
differentially transcribed glycogenes were uploaded to the 
database to explore regulatory network. In the analysis of 
validation in clinical cancer specimens, the housekeeping gene 
jS-Actin was as the reference and the relative expressions value 
were expressed by ACt = Ct target gene- fi-Actin, which means that 
the higher value of ACt of the glycogene, the lower its 
transcription level. 

Results 

The transcriptional profiling of glycogenes 

The transcriptional profiling of the 84 glycogenes was obtained 
by the Human Glycosylation RT 2 Profiler PCR Array and was 



used to produce a hierarchical clustering scheme (Figure. 1A). A 
total of thirty-three glycogenes were differentially transcribed. 
Among them, sixteen glycogenes generated and altered mature N- 
linked glycans, eighteen glycogenes were concerned in the 
synthesis and processing of O-linked glycans, three glycogenes 
acted on glycosphingolipids (Figure 1B-D). Double counting did 
occur among four glycogenes b3gnt3, hexa, st3gall,st8sia3, which 
took part in the synthesis of two kinds of glycans. A total of 
eighteen mRNA, such as mgal5, si3gall, st8sia4, were up- 
regulated in the HCC cell lines with different metastatic potentials 
(MHCC97L, MHCC97H and HCCLM3) compared with the 
non-metastatic HCC cell line Hep3B, while fifteen mRNA like 
clgaltl, gcent4, manlal were down-regulated (Table 2 and 3). 
st3gall could mediate the sialylation of T antigen and the increase 
of st3gall expression had been shown to be one of the major 
mechanisms responsible for the sialylation of T antigen [21]. 
Sialyl-T antigen was tumor-associated carbohydrate antigen 
whose expression was largely increased in some types of cancers 
and associated with poor prognosis [22]. Polysialic acid was a 
carbohydrate composed of a linear homopolymer of oc2, 8-linked 
sialic acid residues and was with a large amount of negative 
charge. The presence of polysialic acid attenuates the adhesive 
property of neural cell adhesion molecule and increases cellular 
motility. st8sia4 was one of the a2, 8-sialyltransferases which could 
add oligosialic and polysialic acid into various sialylated N- 
acetyllactosaminyl oligosaccharides [23-25]. Amounts of pS-1, 6- 
GlcNAc branched N-Glycan were commonly increased in 
malignancies and correlated with disease progression. mgat5 was 
required in the biosynthesis of (5- 1 , 6-GlcNAc -branched N-linked 
glycans attaching to cell surface and secreted glycoproteins. 

Identification in human HCC cell lines and their 
orthotopic xenograft tumors by qRT-PCR 

Seven glycogenes: clgaltl, galnt3, gcnt3, manlal, mgat4a, 
mgal5, mgat5b, were selected to identify the results of the PCR 
Array by qRT-PCR in human HCC cell lines (Hep3B, 
MHCC97L, MHCC97H, HCCLM3) and their orthotopic xeno- 
graft tumors. The housekeeping gene fi-Actin was as the reference 
and normalized the obtained fold changes of differential 
glycogenes. The results of the human HCC cell lines were shown 
in Figure 2A, C, E, G, I, K, M, meanwhile the results from the 
orthotopic xenograft tumors were in Figure 2B, D, F, H, J, L, N. 
clganltl, manlal, mgat4a, galnt3 were all down-regulated 
whereas gcnt3, mgat5, mgat5b were all obviously up-regulated, 
almost all the transcription levels of the glycogenes had significant 
difference. These results supported the transcriptional profiling 
from the PCR Array. 

K-mean clustering and GO analysis 

The thirty-three differentially transcribed glycogenes were 
classified in six clusters according to their trends by X-means 
clustering (Figure 3, Table S2). Among these, only b4gall5, 
galnt7, galntll were up-regulated gradually while man2a2, 
mgall were down-regulated step by step, which implied that 
changes of diverse glycogenes transcription had little related with 
the level of metastatic potentials. GO teams were used to describe 
three attributes of the differentially transcribed glycogenes 
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M N 



ingttUa expression in HCCcell lines mgalja expression in HCC orthotopic xenograft linnors 




Figure 2. Identification of the differentially transcribed glycogenes in human HCC cell lines and their orthotopic xenograft tumors 
by qRT-PCR. A, C, E, G, I, K and M were the HCC cell lines with different metastatic potentials (MHCC97L, MHCC97H, HCCLM3) with the comparison 
to the no-metastatic Hep3B cell line, respectively. B, D, F, H, J, L and N were the orthotopic xenograft tumors of the HCC cell lines MHCC97L, 
MHCC97H, HCCLM3 comparison to the orthotopic xenograft tumor of the Hep3B cell line. The expression of these glycogenes were normalized 
against endogenous mRNA of the housekeeping gene, fl-Actin. 
doi:1 0.1 371 /journal.pone.01 07941 .g002 
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ABC 




Figure 3. K -means clustering of differentially transcribed glycogenes. Each line represented one differentially transcribed glycogene, the 
shape of lines displayed the expression trend of the differentially transcribed glycogenes in the no-metastatic HCC cell line Hep3B and the HCC cell 
lines with different metastatic potentials (MHCC97L, MHCC97H and HCCLM3). 
doi:1 0.1 371 /journal.pone.01 07941 .g003 
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Figure 4. GO analysis of differentially transcribed glycogenes. The differentially transcribed glycogenes were described in 3 categories: 
biological process, subcellular location, molecular function. GO teams were used to describe three as their attributes of the differentially transcribed 
glycogenes. 

doi:1 0.1 371 /journal.pone.01 07941 .g004 
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Table 4. Upstream regulators and target molecules in the regulatory network. 



Upstream Regulator 


EntrezlD 


Molecule Type 


p-value of overlap 


Target molecules in dataset 


ACSL5 


51703 


enzyme 


0.0246 


galnt13 


ADCYAP1 


116 


other 


0.0467 


galnt7, mania! 


AIRE 


326 


transcription regulator 


0.0023 


mgatl, mgat5 


BAK1 


578 


other 


0.0416 


edeml 


CCND1 


595 


other 


0.0418 


galntl 1, st6galnad 


CLOCK 


9575 


transcription regulator 


0.0186 


galntl 1, pomgntl 


ETS1 


2113 


transcription regulator 


0.0352 


b4galt5, mgat5 


ETV4 


2118 


transcription regulator 


0.0416 


b4galt5 


HIST1 HI A 


3024 


other 


0.0003 


b4galt5, man2b1, ogt 


HIST1 HIT 


3010 


other 


0.0003 


b4galt5, man2b1, ogt 


IDH1 


3417 


enzyme 


0.0416 


edeml 


KCNIP3 


30818 


transcription regulator 


0.0232 


edeml 


NCF4 


4689 


enzyme 


0.0146 


st3gall 


SDC1 


6382 


enzyme 


0.0117 


galnt3 


SIRT1 


23411 


transcription regulator 


0.0125 


edeml, mgatl 


STAT5B 


6777 


transcription regulator 


0.0229 


hexa, st8sia4 


TFEB 


7942 


transcription regulator 


0.0445 


hexa 


TMBIM6 


7009 


other 


0.0189 


edeml 


XBP1 


7494 


transcription regulator 


0.0340 


edeml, edem2 



doi:1 0.1 371 /journal.pone.01 07941 .t004 



(Figure 4). The differentially transcribed glycogenes were closely 
related with some types of biological process, such as protein 
amino acid glycosylation, glycoprotein metabolic process, glyco- 
protein biosynthetic process. The subcellular distributions of these 
glycogenes were enriched in golgi apparatus, golgi membrane, and 
so on. They were associated with molecular functions, like 
acetylglucosaminyltransferase activity, mannosidase activity, poly- 
peptide N-acetylgalactosaminyltransferase activity. 

Regulatory network by IPA 

All of the thirty-three differentially transcribed glycogenes were 
uploaded to the database to analyze upstream regulatory events 
and build the possible regulatory network by the IPA (Figure 5). 
The molecule type, p-value of overlap and target molecules of the 
transcription regulators in the regulatory network were described 
and listed in Table 4. We found that a transcription regulator may 
control more than one glycogene, such as SIRT1 which can act on 
edeml and mgatl . At the same time, it was not difficult to detect 
that a glycogene may be simultaneously affected by two or more 
transcription regulators. Instances were described as following: 
galntl 1 was direct regulated by CLOCK, meanwhile it was 
indirectly regulated by CCND1. Moreover, edeml was directly 
influenced by XBP1, KCNIP3 and SIRT1, in the meantime, it 
was indirectly influenced by TMB1M6, IDH1, BAK1. As 
expected, results showed that the different transcription regulators 
tended to have different degree of regulation on glycogene 
expressions and these influences could be measured by the p- 
value of overlap. The smaller p-value of overlap there was, the 
greater influence of the regulation there was. For a given 
glycogene, which of the transcription regulators) had more 
influence could be discriminated by the p-value of overlap. Like 
edeml as a glycogene was influenced by XBP1, SIRT1, TMBIM6, 
IDH1, BAK1, KCNIP3, nevertheless, KCNIP3 and TMBIM6 
play much more important roles. For the seventeen glycogenes in 



the regulatory network, NCF4, SDC1, SIRT1, CLOCK, 
TMBIM6 had smaller p-values and it suggested that they may 
play more important roles in the transcriptional regulation of the 
differentially transcribed glycogenes. 

Validation in clinical cancer specimens by qRT-PCR 

Four differentially transcribed glycogenes, including galnt3, 
gcnt3, manlal, mgatSb were selected to detect the expressions in 
clinical cancer specimens and their results were shown in Figure 6. 
Although the transcription levels of galnt3 and manlal in 
metastatic HCC clinical cancer specimens were lower than in 
the non-metastatic ones, there were no significant differences 
(P = 0.246, P= 0.108, respectively). While the divergent trends of 
gcnl3 and mgat5b between non-metastatic and metastatic HCC 
clinical cancer specimens were statistically significant (P = 0.002, 
P = 0.04, respectively) and they were all up-regulated in metastatic 
HCC clinical cancer specimens, which trends were as same as in 
human HCC cell lines and their orthotopic xenograft tumors. 
gcnt3 (glucosaminyl (N-acetyl) transferase 3, mucin type) is 
required to form GlcNAcp 1 1^6Gal/GalNAc (Figure 6E) [26] 
and mgat5b (mannosyl (a-l,6-)-glycoprotein (3-1,6-N-acetyl-gluco- 
saminyltransferase, isozyme B) acts on the GlcNAc pi,2-Manotl- 
Ser/Thr moiety, forms a 2,6-branched structure in brain O- 
mannosyl glycan (Figure 6F) [27]. 

Discussion 

A growing number of researches indicated that the structural 
modifications of cell membrane glycans had correlation with 
disease progression, metastasis, and poor prognosis [28,29]. 
Essentially, cancer could be regarded as a molecular disease of 
cell membrane glycans because cancer cells were characterized by 
aberrant glycosylation of the surface membrane [30] . Abnormal 
glycosylation was also highly associated with HCC progression. 
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Figure 5. Regulatory network of glycogenes by IPA. The network of glycogenes was derived from the thirty-three differentially transcribed 
glycogenes. The grey ones represented they were included in the differentially transcribed glycogenes whereas the white ones represented they 
were not. There were six types of relationship in the network, A: activation, E: expression, M: modification, PD: protein-DNA interaction, PP: protein- 
protein interaction, T: transcription. Full lines meant a direct action between two nodes, while the dotted lines meant an indirect relationship 
between two nodes. 
doi:1 0.1 371 /journal.pone.01 07941 .g005 



The alterations of glycan branched structures of glycoproteins, 
such as GlcNAc-branched N-glycans, sialic acid and core 
fucosylation, were demonstrated to be highly expressed in HCC 
and these changes were closely related to the transcriptional and 
translational levels of the glycogenes [31-34]. It was generally 
known that AFP-L3, which was the LCA-bound fraction of AFP, 
as a new generation of tumor marker for HCC had been widely 
used in clinic. It had been reported that malignant liver cells 
produce AFP-L3, even when HCC was at its early stage [34]. In 
the meantime, ful8 was known to be a key enzyme of fucosylation 
and it had been implicated in the production of AFP-L3 [35,36]. 
This example suggested that some abnormal glycosylation and the 
differentially transcribed glycogenes might be potential markers of 
biologic malignancy of HCC and maybe they could be used for 
identifying subtypes of HCC and elevating the possibility of 
individual treatment, hence the further research was needed. 



In current study, we obtained and validated the transcriptional 
profiling of glycogenes. This was the first study to select the 
transcriptional profiling of glycogenes associated with HCC 
metastasis by PCR Array and it was a comprehensive analysis of 
glycogenes with relevant N-linked glycosylation, O-Linked glyco- 
sylation, glycosphingolipids in human HCC cell lines, orthotopic 
xenograft tumors and clinical cancer specimens. 

The metastatic potential of HCC cell lines in our experiments 
were strengthening step by step. Hep3B, the non-metastatic cell 
line from ATCC, was established from a patient of 8-year-old 
black boy. Nevertheless, MHCC97L, MHCC97H and HCCLM3, 
all of them can metastasize, whose metastatic potential were 
increasing in general [35] . The original source of these cell lines 
was MHCC97 which derived from the LCI20 tumor line (male 
HCC) and their lung metastatic rate were 40%, 100% and 100%, 
respectively. The node rate of HCCLM3 was 60% [36]. We found 
thirty-three differentially transcribed glycogenes by the Human 
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Figure 6. Validation in non-metastatic and metastatic tissue samples by qRT-PCR. A, B, C, D were the transcription levels of galnt3, gcnt3, 
manlal, mgat5b, respectively. The housekeeping gene, ji-Actin was used to normalize the expression levels in the subsequent quantitative analyses, 
and the higher value of ACt of the glycogene, the lower its transcription level. E, gcnt3 involved in the glycan biosynthetic pathways. F, mgat5b was 
correlated with the pathway of glycan biosynthetic. 
doi:1 0.1 371 /journal.pone.01 07941 .g006 



Glycosylation RT 2 Profiler PCR Array with differences greater 
than 2-fold (PS0.05) as a criteria. A glycomics study using 
glycogene microarray which contained 115 genes by Kang et.al 
reported eighteen glycogenes were up-regulated in high metastatic 
potential HCCLM3 cell lines in comparison to Hep3B cell lines, 
while eleven glycogenes were down-regulated [4]. Guo et.al's 
study using qRT-PCR for quantification found that six glycogenes 
b3galtl, fut8, gal3st2, mgat5, sl3gal5 and sl6gall were overex- 
pressed in MHCC97H cell lines compared with those in 
MHCC97L cell lines; conversely, five glycogenes, b3gni6, 
b3gnl7, galnt6, mgat3 and st8sia2 were down-regulated [7]. 
However, Kang et.al only compared the high metastatic potential 
HCC cell lines HCGLM3 to non-metastasis cell lines Hep3B and 
Guo et.al only compared two metastasis cell lines (MHCC97L and 
MHCC97H). Comparisons both of them carried out were 
between two cell lines, while our results were more confident 
than theirs with four cell lines to be compared, and all of the 
differently transcribed glycogenes had been identificated in 
orthotopic xenograft tumors and validated in clinical cancer 
specimens. Kang et.al used glycogene microarray, by contrast, we 
used the PCR-Array. It was cheaper than by gene microarrays and 
its experimental conditions could be easily controlled [37]. In 
contrast with Guo et.al's result by qRT-PCR, our findings 
obtained more differentially transcribed genes. Additionally, apart 
from the glycogenes which have been reported previously like 
mgall , mgat5 and st3gall, some glycogenes such as "UDP-N- 
acetyl-oc-D-galactosamine: polypeptide N-acetylgalactosaminyl- 
transferase" family members: galnt3, galnt4, galntll, galntl3, 
and "mannosidase, oc class" family members, for instance, 
manlbl, manlcl, man2a2, their different expressions were 
associated with HCC metastasis are the first reported. This 
research revealed that glycogenes which could generate and alter 
sialylated glycan structures or (3-1, 6-GlcNAc branched N-Glycan 
were significantly up-regulated in the metastatic HCC cell lines 
and their orthotopic xenograft tumors. Sialic acids could prevent 
cell-cell interactions through charge repulsion effects [38]. Tumor 
cell lines with more invasive ability expressed more heavily 
sialylated glycan structures. mgat5 could advance tumor growth 
and metastasis by destructing the extracellular matrix or boost the 
production of angiogenic factors by a metal ion-dependent serine 
protease. mgat5 could also promote tumor metastasis by catalyzing 
sugar chain which can change the function of angiogenic factors or 
by directly increasing the gene transcription of angiogenic factors 
[39]. 

The qRT-PCR analysis of the transcribed glycogenes in the 
human HCC cell lines and their orthotopic xenograft tumors 
provided evidences of the same changing trends in statistics (P< 
0.05) to support the analysis of PCR-array results. Nevertheless, 
we found that the respective fold-change value of comparison 
results within metastatic HCC cell lines and their orthotopic 
xenograft tumors had no obvious gradually change trend. It 
revealed that the transcriptional of glycogenes was highly related 
with metastasis no matter whether their metastatic potentials were 
high or low. Therefore, we suggested that our future studies about 
glycosyltransferases should more focus on comparison of whether 
these cells or orthotopic xenograft tumors were able to metastatic 
or not rather than the difference of their metastatic potentials. 



The regulatory network of part of the differentially transcribed 
glycogenes was obtained by the IPA, and it included six types of 
relationship and nineteen upstream regulators. Aside from 
regulators, current data suggests that miRNA is related with the 
expressions of glycogenes. Human and rodents HCC cell lines 
were transfected with miR-122 and miR-34a and the result 
demonstrated that forced expression of these miRNAs was able to 
induce a decrease oifut8 levels and also to affect core fucosylation 
of secreted proteins [40]. The passenger strand miR-17-3p 
repressed expression of galntl , which regulates metabolism of 
liver toxin galactosamine [41]. Further analysis including glyco- 
genes regulation may elucidate the changes of the expressions of 
glycogenes. 

The validation in clinical cancer specimens showed some 
differences of the transcription levels of glycogenes in clinical tissue 
samples were not significant. It might be caused by the part of the 
differences among cell lines, orthotopic xenograft tumors and 
clinical tissue samples themselves. But it was important to note 
there were still a substantial part of glycogenes which were 
differentially transcribed between non-metastatic and metastatic 
tissues. This evidence suggested there might be some potential 
novel drug targets and biological markers to predict metastasis of 
HCC in the differentially transcribed genes. To validate the 
transcription levels of the remaining differentially transcribed 
glycogenes, further investigation is needed and this part of the 
work is undergoing. Moreover, it could not be neglected that our 
study is the first time to provide evidences in human HCC cell 
lines, orthotopic xenograft tumors and clinical tissue samples that 
the up-regulated of gcnt3 and mgat5b were associated with HCC 
metastasis. 

In summary, this study identified thirty-three differentially 
transcribed glycogenes by the Human Glycosylation RT 2 Profiler 
PCR Array. A total of seven glycogenes, including clgaltl , galnt3, 
gcnt3, manlal, mgat4b, mgat5, mgat5b, were further identified by 
qRT-PCR in HCC cell lines and their orthotopic xenograft 
tumors. All of these differentially transcribed glycogenes were 
analyzed by AT-means and GO analysis, and were uploaded to the 
IPA database for regulatory network analysis. What was more, 
galnt3, gcnt3, manlal and mgal5b were selected to further 
validation in clinical cancer specimens and this is the first report 
that the up-regulated of gcnt3 and mgat5b were associated with 
HCC metastasis. The differentially transcribed glycogenes could 
be the potential targets for the novel drugs and might act as the 
potential biological markers to predict metastasis of HCC. 
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